Coarse grained approach for universality classification of discrete models 
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Discrete and continuous models belonging to a universality class share the same linearities and 
(or) nonlinearities. In this work, we propose a new approach that allows for calculating coarse 
grained coefficients and hence to determine the universality class of unrestricted and restricted 
discrete models. Therefor we apply small constant translations in a test space and show how to 
obtain the coarse grained coefficients from the transformed average interface growth velocity. Using 
the examples of the ballistic deposition (BD) model and the restricted solid-on-solid (RSOS) model, 
both belonging to the Kardar-Parisi-Zhang (KPZ) universality class, we demonstrate how to apply 
our approach to calculate analytically the corresponding coefficients of the KPZ equation. Our 
analytical nonlinear coefficients are in agreement with numerical results obtained by Monte Carlo 
simulations. In addition to the BD and the RSOS we study a competitive RSOS model that shows 
crossover regimen between the KPZ and Edwards- Wilkinson universality class. 
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I. INTRODUCTION 

Discrete and continuous models belonging to the same 
universality class are usually characterized by a unique 
set of power-law exponents related to each other by scal- 
ing relations [H-I^. Beyond that, discrete and continuous 
models within the same class share the same linearities 
and (or) nonlinearities [5j]. All this properties are result 
of the evolution rules determined by the transition rates 
between surface configurations. Usually this rules are 
based on nested "if-then-else" structures, which alterna- 
tively can be expressed by means of Heaviside unit-step 
functions 8. Starting from the Master equation of the 
time evolution of the configuration probability, it is pos- 
sible to derive an associated discrete Langevin equation 
performing a Kramers- Moyal expansion Please 
note that the Q -function is a distribution or generalized 
function and therefore the drift function of the discrete 
Langevin equation is a distribution, too. The continuous 
Langevin equation can be achieved from its discrete coun- 
terpart following the approach introduced by Vvedensky 
et al. [§[ . This approach is based on regularization tech- 
niques and on the coarse grain approximation, consistent 
with the limit of lattice constant tending to zero 
The regularization method involves the replacement of 
the Heavisides Q by smooth functions 9^. Here the reg- 
ularizing parameter e has to be chosen in a way that 
9^ ^ Q when e — s> 0+. The coarse grain approximation 
is carried out by a Taylor expansion around the origin of 
either the entire regularized drift function or each single 
regularized function 9^ and its arguments (e.g. Lapla- 
cian or gradient). It can be observed, that the results 
of this two expansion methods are ambiguous, since they 
strongly depend on the chosen regularization [loj . 
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Within the group of discrete models, which describe 
deposition or evaporation processes, the most studied 
ones belong to the Edwards- Wilkinson (EW) or the 
Kardar-Parisi-Zhang (KPZ) universality class. Models of 
the EW class treat slow deposition processes with relax- 
ation under gravity Models of the KPZ class involve 
deposition (or evaporation) processes that permit lateral 
growth (or decrease) [ll] . Usually studies of these models 
were focused on the characterization of statistical surface 
growing properties, e.g. roughness and growth exponents 
of the surface width and scaling function in the steady 
state. The KPZ universality class includes several dis- 
crete models, e.g. the restricted solid-on-sohd (RSOS) 
model [13], the Eden model [H, [l^, ballistic deposi- 
tion (BD) model [H [23| , and the deposition-evaporation 
model [2l[. Among these models, many authors chose 
to study the RSOS model, as it shows a very good scal- 
ing behaviour even with small systems and in addition 
reaches rapidly the steady state. First, Park and Kahng 
[H have derived the KPZ equation [3| from the RSOS 
model. More recently using the same method, Oliveira 
et al. [lH has derived the KPZ equation from a com- 
petitive restricted solid-on-solid (CRSOS) model involv- 
ing deposition and evaporation, with probabilities p and 
1 — p, respectively. Other authors have derived the KPZ 
equation starting from the BD model following different 
schemes of regularizations [13, ■ 

In this work we propose a new approach that allows 
for deriving the linear and (or) nonlinear coefficients of 
a corresponding continuous equation starting from a dis- 
crete model, and thus to classify discrete models. We 
introduce the concept of interface configuration space 
(ICS) in which discrete processes are defined. The com- 
ponents of a point in the ICS are the height differences 
between neighbour columns and the site column evolved. 
Within this configuration space the drift function of dis- 
crete Langevin equation and the stationary probability 
density function (spdf) are defined. Additionally, we as- 
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sume that the spdf for a restricted (or unrestricted) dis- 
crete process is embedded in a space of test functions 
with compact support (or rapidly decreasing behaviour). 
We show that proper constant translations within the 
test space applied to the test function transform the av- 
erage velocity of the interface growth. The coefficients 
that occur under the proposed transformations corre- 
spond to the coarse grained coefficients of continuous dif- 
ferential equation which defines the universality class. In 
the framework we focus on the RSOS and the BD model 
of the KPZ class, since they can be taken as a reference of 
a restricted and unrestricted processes respectively where 
our theoretical approach can be applied. 



(2) 

Here Qij = K^^' is the noise intensity. Introducing the 
probability current [dj = d/dhj) 

J, = K,P-\dj{Q,jP) , 

the Fokker-Planck equation ([T]) is 

dtP = -a, J, . 

In the steady state dtPst = 0, where P^t = Pst{<Zj) is the 
spdf with the property J2(rez'f Pstiz) = 1 • 



B. Transformations on interface configuration 



A. Prom Master equation to discrete Langevin 
equation 

Let us consider a surface configuration H, which is 
determined by a set of heights {hj} corresponding to 
the columns j. The transition rate W{H,H') between 
two surface configurations H and H' , for a process that 
evolves (increasing or decreasing r^-units) at the selected 
columns fc, is 



1 



N 



W{H, H') = -Y. ""^(S, H')Aih'„ hk+Tu) n Hh'j.h,) , 



k=l 



where A(a;, y) is equal to 1 if a; = y and equal to oth- 
erwise. Here the hopping rules ujk{H,H') take the re- 
straints process into account. For restricted models, e.g. 
the RSOS U!k is function of the height difference between 
the selected column k and its next neighbours. For un- 
restricted models, e.g. the BD, cj^ = 1. The growth 
rk of the columns k depend on the underlying discrete 
growth model, e.g., Vk = 1 for the RSOS model and 
rk = max(fTfc^,crfc__,l), with ak± = hk±i - hk for the 
BD model. The first and second transition moments are 



The first approach that allowed for characterizing non- 
linearities of discrete models was the surface tilt trans- 
formation [H, [IJ . It is based on a global transformation 
V/i — )■ V/i -I- s which is applied to continuous KPZ equa- 
tion [h = /i(x, t) with X e M"] 

dh 

'dt ' • ■ 2 
and transforms the average velocity of interface growth 
asi 



^ |V/ip +77(x,t) , 



(3) 



i>(0) -J- w(s) = u(0) + 



(4) 



with s = |s|. Employing the tilt transformation and 
Monte Carlo techniques, several authors characterized 
nonlinearities of discrete models belonging to the KPZ 
universality class [23|-[25|. In this work we propose di- 
rect transformations of the discrete Langevin equation 
([^. For simplicity we consider only two next neigh- 
bours (Tj — {uj^^cTj ) and use the tilt transformation 

hj — > hj -I- sj with s g M. In this case the hight dif- 
ferences between the next neighbours are transformed as 



<7j^ ± s. We define the two auxiliary variables 



H' 
H' 



2 \'^]+ ^i- 1 

'j = \ W3+ + CTj-) 



\{hj+i 
\{hj+i 



hj-i) , 
2hj + /ij-i) 



{K-h,)ih'-h,)W{H,H') 



This new variables are the standard (or post-point) dis- 
cretization of the gradient 13[ and the discretization of 
the Laplacian. Using the tilt transformation, they trans- 
where aj — {a^i, . . . , ajN) G and ajk — hk~ hj is the form as rjj ^ r/j + s and Cj ^ Cj i-e- the transfor- 

height difference between the neighbour column k and 
the selected column j. Using a Kramers- Moyal expansion 
the following Fokker-Planck equation is derived from the 
Master equation [P = P({(Tj},t)] 



dP 
'dt 



2 dh.dhj^ 



{K\fP) . (1) 



The corresponding Langevin equation is {dt = d/dt) 

dth,=K,{<j,)+rj,{t) , (2) 



mation shifts the "gradients" and leaves the "Laplacian" 
invariant, respectively. In contrast, the transformation 
aj^ (Tj^ + r leaves "gradient" invariant and changes 
the "Laplacian" to Q (j + r. These transformations 
may be interpreted as translations along the "gradient" 
and "Laplacian" directions of the surface configuration 
space. 



II. DISCRETE PROCESSES REVISITED FROM 
THE THEORY OF DISTRIBUTIONS 



is Gaussian 



where Kj = K^^^ is the drift term and rjj 
white noise with {rjj) = and {r]i{t)rij{t')) = Qij S{t — t') 



Let Z^= {ae Z'^/Pstiz) ^ 0} be the set of surface 
configurations in the steady state. Here is the lat- 
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tice in TZ C consisting of points with integer coordi- 
nates. The drift term of the discrete Langevin eq. ^ can 
be assumed as a generahzed function K(q) with q G 

evaluated at g = (Tj <E . In order to calculate statisti- 
cal observables of discrete processes, with and without 
restrictions, it is needed to define test functions ip on 
which the generalized functions of g are applied. For re- 
stricted processes applies: ip G 'D{TZ^). Here V is the test 
space of C°°-functions with compact support in TZ^ = 
e - neighborhood(7^^) = e R^/d(£,7^^) < e} , 

where C (K^ \ Z^) U is a compact set. The 
support of the function ip : TZ^ — !■ C, denoted supp(93), 
is the closure of the set {g e TZ^ /fig) 0} • The com- 
pact support excludes interface configurations in Z^\Z^ 
since supp(</3) C TZ^ . In contrast, unrestricted processes 
require that (p G S{R^). Here S is the test space of C°°- 
functions that decay and have derivatives of all orders 
that vanish faster than any power of g~^ (a = 1, . . . , N). 
Our notation about distributions is the standard, first 
established by Schwartz [2^ . The application of a distri- 
bution f € V (dual space of 2?) to test function (p E V 
is defined by 

/ fig) fig) d^'g. (5) 

Please note that (/, ip) is the expectation value of / us- 
ing the test function ip as real analytical "representation" 
of the spdf Pst- Thus, the test function is normed, i.e. 
(1,</?) = 1 . The translation Tq, of a distribution /, de- 
noted Taf or simply fa, extends the definition given by 
eq.® 

{Taf,f) = {f,T^af), (6) 

where the translation operator is defined by 
T^: y t-^ y — xify ,x d |27| . We assume that the test 
function ip takes fixed values in the configuration space 
given by the spdf P^t, i.e. {S^ , ip) = ipia) = Pst (2) 
for all a G , where — T„5 and 5 is the Dirac 
distribution. Applying a translation T„ to a point 
g G supp((p), the test function transforms as — )■ T^ip if 
ig~1i) G supp((/5) and the average velocity of the growth 
interface change as u(0) — ?> with 

viu) = {K , Tuf) = / Kig) ipig - u) d^g . (7) 

For small translations the Taylor expansion of ipig — u) 
around w = is 

fig-u) = <p>ig)-da'P>\^^^ ^c, + l9lMu=0 UaUj3 + --- . 

(8) 

Here repeated subscripts imply sums, Ua is the a-th com- 
ponent of u, da^d/dga, and d'^n ^ d"^ / idgadgp), with 



a = j3 = 1, . . . ,iV. Since the test function ip is known 
only at points of the lattice Z^ , its derivatives can not 
be calculated explicitly. In contrast, the drift distribu- 
tion K is derivable in all points. Since the test function 
either has compact support or decreases rapidly one can 
take advantage of the following identity 

(x,5v...^)-(-irfe...x,^). (9) 

Using eqs. dlHS]), the average velocity is transformed ac- 
cording to 

v{ll) = (^K, ip'j + (daK, ip^ Ua + \(dafiK, ip^ UaUp H . 

(10) 

The physical meaning of the translations in test space can 
be made clear considering the case N = 2. It corresponds 
to a discrete process that describes the evolution of an 
interface in (l-l-l)-dimensions. -ft'(f?i,£'2) is symmetric 
under the exchange gi -H> g2, for an isotropic process. 
Consequently (pigi, P2) is symmetric under this exchange, 
too. Considering a translation u in the supp((/3) with 
components ui — r + s and U2 — r — s, the lower orders 
in the Taylor series expansion are 

vir,s)^(^K,ip) + (^idi+d2)K,ip)r 

+ ^(^id!,-2d!, + di,)K,ip) s\ (11) 

Here the linear term in s vanishes because (9i — d2)K is 
antisymmetric. Taking into account that ((3i +d2)K and 
idfi — 2 + 822) K are symmetric eq. ([TT|) is 

vir,s)~VQ + iyr+^Xs'^ , (12) 

with the coefficients (a, /3 = 1, 2 with a ^ /3) 

1/ = 2 , (^) , (13) 

X^2(^{dla~dlf,)K,ip) , (14) 

vo^(K,p>)^viO,0) , (15) 

where repeated subscripts do not imply sums. Assum- 
ing that higher-order nonlincarities can be neglected, dis- 
crete models belong to KPZ universality class (A ^ 0) if 
(^Qa ~ ^a;?)^ "^^^^ symmetric under the ex- 

change Pi ^ ^2 • Otherwise, the discrete models belong 
to Edwards- Wilkinson (EW) universality class (A = 0). 
Equation (|12p can alternatively be obtained applying 
translations to the gradient and the Laplacian term of 
the KPZ equation, i.e. V^/i V'^h+r and Vh Vh+s, 
and then taking averages over the interface. The set of 
parameters {v,X} defined by eqs. (fT3|) and ((T4)) is unique 
for each isotropic discrete model. 

III. RSOS MODEL 

For a single-step aggregation process in one dimen- 
sional substratum, the RSOS condition is \gk\ < 1 with 
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k = 1,2. Evolving from a flat interface we observe 
9 interface configurations defined by = G 
'^'^ / i;i = ~ljO, 1}. The test function 1^9(^31, Q2) is taken 
as a C°°-function with compact support, defined in TZ^ = 
e — neighborhood(7?.^) an open subset of M? that ex- 
cludes all points of 1? \ Z'^ . For simplicity we take 
7^2 = G RVIftI < 1, = 1,2} (see fig. ^. The 
drift of Langevin eq. ^ is the generalized function 



(16) 



where we take a/r = 1 without loss of generality and the 
Heaviside function Q{z) is equal to 1 if z > and equal 
to otherwise. The derivatives of the K up to 2-nd order 

are (a, (i — 1,2 and a ^ (3) 



(17) 



Applying these distributions to tp we obtain the coeffi- 
cients V and A defined by eqs. ([13]) and (IT4)) . respectively. 
Explicitly 



v = 2 I (^(0,y)dy, 







¥'(0,0) 



(0,y) 



dy 



(18) 
(19) 



Numerical results: In order to evaluate eqs. (|18p and 
p9)) we taking e ^> 0+ and assume that (/s and 9a;(/7 at 
X — Q are very smooth functions oi y in (0, 1). Straight- 
forward we obtain v ~ </'(0,0) + (^(0,1). We obtain 
A ~ -2 ((^ 9a;<y9 + i 9^j^V')J(o,o) replacing d^ip\i^o^y) ~ 
^2 . /[Y9|_ Please note that 



c'2,<pJ(o,o) + 2/c)4<pJ(o,o) in eq. 
5a;'/'J(o.o) - [¥'(1,0) - (^(-l,0)]/2, and 9:cy¥'J(o.o) - 
-[(^(1,1) - 2v3(-l,l) + (y3(-l,-l)]/4. With the val- 
ues given in the caption of fig. [1] we obtain v ~ 0.2102 
A ~ -0.2171. 



IV. CRSOS MODEL 

The simplest competitive model that shows a crossover 
between EW and KPZ universality classes is the 
deposition-evaporation model with RSOS condition. The 
deposition (or evaporation) happens with probability p 
(or 1 — p), in a one dimensional substratum. This model 
was first numerically studied by Amar and Family far 
from crossover [25]. Later Olivera et al. showed analyt- 
ically the X-p relation at the associated KPZ equation 
[11] . Both, the interface configuration space Z^ and test 
space TZ"^ are identically to the ones of the RSOS model. 
The drift term of the CRSOS model is 



Kigi, Q2) ^pQ{Qi) Q{Q2) - (1 - P) Q[-Qi) Q[-Q2) 




FIG. 1. For the RSOS model we plot a smooth test 
function as function of {qi,Q2)- The non zero values 
of the spdf (bullets) obtained by Monte Carlo simulation 
are: ^3(0,0) = ^5(1,1) = 0.1158, ^5(0,1) = (/^(l.O) = 
0.0944, <^(-l,l) = V5(l,-1) = 0.1249, <^("1,0) = 
V5(0, -1) = 0.1198, and (^(-1,-1) = 0.0903. See that 
supp(<^) is the projection of <j9 7^ on the plane (light gray). 



and its derivatives up to 2-nd order are 

d^K = 5{q^)[pQ{q0) + {I - p)Q{-gp)\ , 
d^pK = {2p-1)5{qo,)5{q,3), (21) 
d^o^K = 5'{q^) [pQ{qp) + (1 - p) e(-e^)] . 

Applying these distributions to if we obtain the KPZ 
coefficients defined by eqs. (|13p and HH). The nonlinear 
coefficient is [281 



A = -(2p- 1 



2</'(0,0) 



sgn(y) d^ip 



dy 



dx^p]^ dy 
i-e -I (o,y) 



(22) 



Thus, when p = 1/2 the process belongs to the EW uni- 
versality class coinciding with the results obtained by 
other approaches [Til . [29| . Please note that A = iff 
dx'p\{o,y) is odd function of y, which happens only for 
p — 1/2 as observed by symmetry of the process. To 
understand the latter just note that K in eq. (|20p with 
p = 1/2 is invariant under the exchange gj —gj with 
j = 1,2. When this symmetry breaks the KPZ univer- 
sality class appears. 



V. BALLISTIC DEPOSITION MODEL 

The evolution of the ballistic deposition (BD) model in 
(l+l)-dimcnsions is defined for all (^1,^2) € by the 
distribution 



K{gi, g2) = max(ei, ^2, 1) , 



(23) 



where we define a/r = 1 without loss of generality. Al- 
ternatively the BD evolution rules can be written as 



(20) 



if = 56(5-1) + [1-6(5-1)] , 



(24) 



FIG. 2. For the BD model we plot the drift distribution K 
as function of {qi-,Q2), given by eq. (I23p . The bullets show 
the discrete drift K{i,j) where € 1? ■ Dashed line shows 
where the distribution is not smooth. 



FIG. 3. For the BD model we plot a decreasing test function ip 
as function of (gi, £12), which take fixed values (pii,j) obtained 
from the Monte Carlo simulation (bullets). At the origin is 
smooth though not visually fit. 



with 



9 = max(gii,£<2) = C+\v\ 



(25) 



Here the auxiliary variables C and 77 are again defined as 
C = '^{91 + 92} and 7] = ■^{91 — 92)- The 2-nd order Taylor 
expansion of K{(^ + r,ri + s) around r = and s = is 

K{C + r, ry + s) ~ KiC, rf) + [rdc^ + sd,^)K 

+ l{r^dl^-2rsdl^ + s^dl^)K . (26) 

Taking into account that ~ \ri\' — sgn(?7) = sgn(£)i — 
£12) [13 and d'^^g = 2 5{vi) ^ A5{qi- Q2), where sgn(?7) = 
0(77) — 6(— 77) is the sign function, one achieves 



d^K = {g - 1) 5{g - 1) + Q{g - 1) , 

(.9 -1)5(5-1) + 9(5-1)] sgn(77), 

(5-1)% -1)1 +<5(5-l), 



dr,K = 



dt,^K = 2 



(5-l)5(5-l)J +<5(5-l) sgn(77), (27) 
(5-l)%-l) + e(5-l)l <5(77) 



(5-l)%-l) sgn(r;) 



+ %-l)[sgn(r;)]2. 

Applying the distribution [eqs. (|26p ] to the test function 
(f one obtains eq. (fT2|) with 

v^id^K,^) - (e(5-i),(/^ 

\={dl^K , v^) - 4 (e(5 - 1) 5{q^ - ^2) , 

+ (5{g-l)[\-A{Q^,g2)],ip), 
VQ^{K,^) = l-v+{gQ{g~l),ip). (28) 



Here we state only non zero contributions and 
^{qitQ2) — 1 — [sgn(gi — ^2)]^- Note that the application 
of the distribution (5 — 1) 5{g — 1) (or its derivatives) on 
a test function is always zero. Also, the antisymmetric 



sign distribution applied to a symmetric test function is 
zero. The term of 2-nd order in r has non-zero coeffi- 
cient (5[g — 1),(^), which is negligible compared to u. 
In order to calculate the coefficients and to understand 
their meaning it is convenient to express the eqs. 
explicitly, e.g. 



A = 4 



+ C30 



^{x, x) da; + 2 



ip{\,x) dx . (29) 



Here e > is a small constant that avoids double evalu- 
ation of the integral at x = 1. One observes, that on one 
hand the coefficient A depends on integrals of in points 
where the distribution K is not smooth (see fig. [5]). On 
the other hand, the coefficient v depends on the integral 
of the if where K is not constant. Numerical simula- 
tions confirm these observations. In fact only the con- 
figurations for all i < 1 and (l,i) or («, 1) for all 
i < \ contribute to A, whereas other configurations of 
lateral growth are negligible. These conclusions were not 
reached by other approaches and could be confirmed by 
other unrestricted models with similar features, i.e. dis- 
crete Langevin equation with not smooth drift functions. 

Numerical results: The values of spdf obtained at 
points of configuration space of rules by Monte Carlo sim- 
ulation in steady-state regime has allowed us to calculate 
the coefficients. The numerical values are A = 0.3265 and 
vq = 2.1356 which are in excellent agreement with those 
reported by Barabasi and Stanley [4| using Monte Carlo 
simulation with initially tilted interface. Additionally, we 
calculated the value v = 1.0302 which is not reported by 
other authors. 



VI. CONCLUSIONS 

In this letter we have proposed a new approach that al- 
lows for calculating coarse grained coefficients and hence 
to determine the universality class of unrestricted and re- 
stricted discrete models. We present a theoretical frame- 
work which is based on well-established concepts of the 
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distribution theory. 

First we introduce the concept of a interface configu- 
ration space (ICS) where discrete processes occur. Ex- 
tending the ICS to the real we define an analytical test 
space. The stationary probability density function is em- 
bedded in a test function space which is generated by test 
functions with compact support or rapidly decreasing be- 
haviour. The drift term of the discrete Langevin equation 
is understood as generalized function or distribution. We 
apply small constant translations in the test space and 
show how to obtain the coarse grained coefficients from 
the transformed average interface growth velocity. The 
coefficients are determined by linear combinations of n-th 
order partial derivatives of a transformed drift distribu- 
tion applied to the test function. 

Then we study exemplarily deposition models without 
relaxation, with and without restrictions, the RSOS and 
the BD models, that define spaces of test functions that 
are rapidly decreasing or with compact support, respec- 
tively. For these (l-l-l)-dimensional models, we deter- 
mine analytically the coarse grained coefficients of KPZ 
class. The numerical values that we have determined for 



the nonlinear coefficient A are in very good agreement 
with those calculated by Monte Carlo simulation using 
the interface tilting method. In addition, our approach 
allows us to calculate the linear coefficient v for these 
two models, which were not known until now. Apart 
from the pure deposition models we also study a compet- 
itive RSOS model that describes a restricted deposition- 
evaporation process. By this example we show how the 
symmetries of the distribution K explain the crossover 
behaviour from KPZ universality class to EW ones. We 
observe that in all study cases A is expressed through in- 
tegrals of the test function or its derivatives on the set 
of points where the K distribution is not smooth or has 
edges or corners. 



ACKNOWLEDGMENTS 

We thank Sandra Molina (Mathematics Department 
at FCEyN-UNMdP) for her useful comments on Distri- 
butions Theory. D.H. thanks CONICET for its support. 



[1] F. Family and T. Vicsek, J. Phys. A: Math. Gen. 18, L75 
(1985). 

[2] F. Family, J. Phys. A: Math. Gen. 19, L441 (1986). 

[3] F. Family, in Springer Proceedings in Physics, Vol. 38, 
edited by R. JuUien, L. Peliti, R. Rammal, and N. Boc- 
cara (Springer- Verlag, Berlin, 1988) p. 193. 

[4] A.-L. Barabasi and H. E. Stanley, Fractal concepts in sur- 
face growth, 1st ed. (Cambridge University Press, Cam- 
bridge UK, 1995). 

[5] Z.-W. Lai and S. Das Sarma, ]Piys. Rev. Lett"] 66, 2348 
(1991). 

[6] H. A. Kramers, Physica 7, 284 (1940); J. E. Moyal, J. R. 
Stat. Soc. B 11, 150 (1949); R. F. Fox, Phys. Rep. 48, 
179 (1978). 

[7] N. G. Van Kampen, Stochastic Processes in Physics and 
Chemistry (North Holland, Amsterdam, 1981). 

[8] D. D. Vvedensky, A. Zangwill, C. N. Luse, and M. R. 
Wilby, Phys. Rev. E 48, 852 (1993). 

[9] K. Park and B. Kahng, Phys. Rev. E 51, 796 (1995). 
[10] M. Predota and M. Kotrla, Phys. Rev. Lett. 54, 3933 
(1996). 

[11] T. J. Oliveira, K. Dechoum, J. A. Redinz, and F. D. A. 

Aarao Reis, |Phys. Rev. E| 7 4, 011604 (2006). 

[12] G. Costanza, 'P hys. Rev. E| 55, 6501 (1997); [Physica A| 

388, 2600 (2009). 
[13] R. C. Buceta, Phys. Rev. E 72, 017701 (2005). 
[14] D. Muraca, L. A. Braunstein, and R. C. Buceta, 

[Ph^. Rev. E 69, 065103 (2004). 

[15] S. F. Edwards and D. Wilkinson. i Proc. R. Soc. Lond. Al 

381, 17 (1982). 



[16] M. Kardar, G. Parisi, and Y. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 

[17] J. M. Kim and J. M. Kosterlitz, Phys. Rev. Lett. 62, 
2289 (1989). 

[18] M. Eden, in Proceedings of the Fourth Berkeley Sympo- 
sium on Mathematical Statistical and Probability, Vol. IV, 
edited by J. Neyman (University of California Press, 
Berkeley CA, 1961) pp. 223-239. 

[19] M. Plischke and Z. Racz, Phys. Rev. Lett. 53, 415 (1984). 

[20] M. J. Void, J. Colloid Sci. 18, 684 (1963); D. N. Suther- 
land, 22, 300 (1966); P. Meakin, P. Ramanial, L. M. 
Sander, and R. C. Ball, |Phys R ctTA 34, 5091 (1986). 

[21] M. Plischke, Z. Racz, and D. Liu, Phys. Rev. B 35, 3485 
(1987). 

[22] J. Krug, J. Phys. A: Math. Gen. 22, L769 (1989). 
[23] D. A. Huse, J. E. Amar, and F. Family, Phys. Rev. A 41, 
7075 (1990). 

[24] J. Krug and H. Spohn, Phy s. Rev. Lett . 64, 2332 (2005). 

[25] J. G. Amar and F. Family, Phys. Rev. A 45, 5378 (1992). 

[26] L. Schwartz, Theorie des distributions, nouvelle ed. (Her- 
mann, Paris, 1966). 

[27] J. F. Colombeau, New generalized functions and mul- 
tiplication of distributions (North Holland, Amsterdam, 
1984). 

[28] Usepe(y) + (l-p)e(-y) = i(2p - 1) sgn(y) + i[l + 
A(j/,0)]. 

[29] A. Lazarides, |Phys. Rev. E|73, 41605 (2006). 

[30] R. P. Kanwal, Generalized Functions: Theory and tech- 
nique (Academic Press, New York, 1983). 



